## # A tibble: 277 × 4
## folder plot_id simul simulation_id
## <chr> <chr> <chr> <int>
## 1 output-cmd_1.txt RETZ_00972_01 simulation_10 1
## 2 output-cmd_1.txt RETZ_00972_01 simulation_11 2
## 3 output-cmd_1.txt RETZ_00972_01 simulation_12 3
## 4 output-cmd_1.txt RETZ_00972_01 simulation_13 4
## 5 output-cmd_1.txt RETZ_00972_01 simulation_14 5
## 6 output-cmd_1.txt RETZ_00972_01 simulation_15 6
## 7 output-cmd_1.txt RETZ_00972_01 simulation_16 7
## 8 output-cmd_1.txt RETZ_00972_01 simulation_17 8
## 9 output-cmd_1.txt RETZ_00972_01 simulation_18 9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19 10
## # ℹ 267 more rows
## # A tibble: 278 × 4
## folder plot_id simul simulation_id
## <chr> <chr> <chr> <int>
## 1 output-cmd_1.txt RETZ_00972_01 simulation_10 1
## 2 output-cmd_1.txt RETZ_00972_01 simulation_11 2
## 3 output-cmd_1.txt RETZ_00972_01 simulation_12 3
## 4 output-cmd_1.txt RETZ_00972_01 simulation_13 4
## 5 output-cmd_1.txt RETZ_00972_01 simulation_14 5
## 6 output-cmd_1.txt RETZ_00972_01 simulation_15 6
## 7 output-cmd_1.txt RETZ_00972_01 simulation_16 7
## 8 output-cmd_1.txt RETZ_00972_01 simulation_17 8
## 9 output-cmd_1.txt RETZ_00972_01 simulation_18 9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19 10
## # ℹ 268 more rows
complete_scene <- import_output_scene( output_name = “complete”, base_path = “C:/Capsis4/data/forceps/clementine/Test_protocole/”, output_file = “output/protocole_completeScene.RData”) %>% filter(date <= “2100”) %>% group_by(plot_id, date, scenario) %>% summarise( totalVolumem3.ha = sum(totalVolumem3.ha, na.rm = TRUE), totalBiomass_t_ha = sum(totalBiomass_t_ha, na.rm = TRUE), .groups = “drop” )
| Plot ID | Surface | Type de peuplement | Essence déterminante | Essence dominante 1 | Essence dominante 2 | Structure et occupation du sol | Surface terrière (m2/ha) | Diamètre moyen (cm) | Âge médian (ans) |
|---|---|---|---|---|---|---|---|---|---|
| RETZ_00230_03 | 2.50 | FHETE | HET | HET | F | 0.0 | 7.5 | 12.0 | |
| RETZ_00816_02 | 0.94 | FHET1 | HET | HET | F | 17.5 | 12.5 | 12.0 | |
| RETZ_00901_01 | 10.23 | FHETG | HET | HET | CHP | F | 21.0 | 57.5 | 144.5 |
| RETZ_00972_01 | 4.59 | FHETG | HET | HET | CHS | F | 31.0 | 57.5 | 144.5 |
| RETZ_01008_05 | 3.38 | FHETS | HET | HET | CHS | F | 0.0 | 3.5 | 12.0 |
| RETZ_01067_02 | 1.73 | FCHPE | CHP | CHP | F | 0.0 | 7.5 | 12.0 | |
| RETZ_01210_04 | 3.10 | FCHSE | CHS | CHS | F | 0.0 | 7.5 | 12.0 | |
| RETZ_01307_02 | 1.28 | FDOUG- | DOU | DOU | F | 29.0 | 52.5 | 59.5 | |
| RETZ_01368_01 | 7.91 | FHETM | HET | HET | CHA | F | 31.0 | 37.5 | 114.5 |
| RETZ_01415_02 | 2.23 | FCHPG+ | CHP | CHP | CHA | F | 17.5 | 62.5 | NA |
Pour rappel des itinéraires :
Je ne suis pas sûre de pourquoi il y a des chutes de surface terrière dans la simulaion 3 (sans coupe), l’itinéraire étant bien : 150_3_0.5_0%_FSyl-80. peut-être juste mortalité naturelle
Sur les 80 ans de simulations on calcul les indicateurs suivants :
Indices de diversité de Hill (richesse, Shannon, Simpson, Hill 0, Hill 1, Hill 2)
Volume extrait par scénario
Biomasse sur pied totale et par espèce